Nature of an intermediate non-Fermi liquid state in Ge-substituted YbRh 2 Si2: 
Fermionized skyrmions, Lifshitz transition, Skyrmion liquid, and Gruneisen ratio 
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We propose a skyrmion liquid state for the non-Fermi liquid (NFL) phase in Ge-substituted 
YbRli2Si2, where skyrmions form their Fermi surface, argued to result from the strongly coupled 
nature between skyrmions and itinerant electrons. The fermionized skyrmion theory identifies the 
antiferromagnetic (AF) transition with the Lifshitz transition, where the quantum critical point 
(QCP) is characterized by the dynamical critical exponent z = 2. Nonlocal interactions between 
skyrmions allow a critical line above the AF QCP, which originates from the Kondo-coupling effect 
with itinerant electrons. This critical line is described by the skyrmion liquid state, which results 
in Landau damping for spin fluctuations, thus characterized by z = 3. As a result, the Gruneisen 
ratio is predicted to change from ~ T -1 at the AF QCP to ~ T~ 2/3 in the NFL phase. 

PACS numbers: 71.10.Hf, 71.10.-w, 71.27.+a 
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I. INTRODUCTION 

Recent experiments on Yb-based heavy fermion sys- 
tems of YbRh 2 (Sio.9 5 Geo.o5)2 0, Yb(Rh . 9 4lr .06)2Sil2 
0, YbAgGe 0, and /3-YbAlB 4 have uncovered 
the appearance of non- Fermi liquid physics over a finite 
zero— T region of the magnetic field- or pressure-tuned 
phase diagram, rather than at a single quantum criti- 
cal point (QCP). The implication of these experiments 
casts a doubt on our fundamental understanding for cor- 
related electrons in metals because the present theoretical 
framework cannot access such a non-Fermi liquid phase 
away from quantum criticality above one dimension in 
the clean system The Fermi liquid phase has been 
regarded as the only fully symmetric quantum ground 
state for metals away from quantum criticality above one 
dimension in the absence of randomness, sometimes re- 
ferred as the mother state of metals @. A non- Fermi 
liquid phase away from quantum criticality, proposed 
to be another stable fixed point in the renormalization 
group sense, is extremely difficult to reach in two or 
three dimensions except for special cases such as the frac- 
tional quantum Hall phase because accessible stable fixed 
points are non-interacting like the Fermi liquid state in 
our perturbative many-body theoretical framework, gen- 
erally speaking. On the other hand, the non-Fermi liquid 
at a QCP is identified with an unstable (either interact- 
ing or non-interacting) fixed point, the origin of which 
is the presence of critical fluctuations associated with a 
quantum phase transition [7[. 

In this paper we revisit an antiferromagnetic quan- 
tum phase transition in the heavy-fermion system, where 
an antiferromagnetic order from localized magnetic mo- 
ments becomes destroyed by the Kondo effect from itin- 
erant electrons. In particular, we focus on the non-Fermi 
liquid phase in Ge-substituted YbRli2Si2, where tuning 
magnetic fields allows such an intermediate state between 
antiferromagnetic and heavy fermion Fermi-liquid phases 



First, we notice the low antiferromagnetic transi- 
tion temperature at zero magnetic field. If we interpret 
the antiferromagnetic transition with the condensation 
of skyrmion excitations as the superconducting transi- 
tion can be understood by the condensation of vortices 
in the dual picture, the low Neel temperature implies 
that skyrmions can be easily created. Then, our problem 
can be mapped into skyrmion dynamics interacting with 
itinerant electrons in the presence of the Kondo-coupling 
effect. This view point reminds us of the phase diagram 
for type-II superconductors in the plane of magnetic field 
and temperature, where rich vortex phases emerge to 
cause non- Fermi liquid physics [8[ . Actually, the analogy 
between our skyrmion-dynamics problem and the vortex 
phase diagram serves an essential physical picture in this 
study. A novel vortex phase was proposed in the region 
above the vortex lattice state, identified with a vortex 
liquid phase, where the strongly coupled dynamics be- 
tween vortices and itinerant electrons gives rise to such 
an exotic liquid state @. 

We propose a skyrmion liquid state for the nature of 
the non- Fermi liquid phase in Ge-substituted YbRli2Si2, 
the characteristic feature of which is the power-law be- 
havior in the skyrmion density-density correlation func- 
tion. A key ingredient in this theoretical proposal is 
to fermionize skyrmion excitations at the antiferromag- 
netic QCP. Although this fermionization procedure can- 
not be justified above one dimension in a rigorous sense 
a priori, essentially the same idea has been applied to 
describe non-Fermi liquid transport phenomena in the 
vortex liquid phase [9J, where the Chern-Simons term, 
which appears from the payment for statistical trans- 
mutation, becomes irrelevant at criticality [10j. This 
fcrmionic skyrmion conjecture identifies the antiferro- 
magnetic quantum phase transition with the Lifshitz 
transition for skyrmion excitations, where the skyrmion 
chemical potential touches the skyrmion band, and thus, 
the dynamical critical exponent is given by z = 2. If 
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we translate tuning of the magnetic field with control- 
ling of the skyrmion chemical potential, we are driven 
to conclude that the Fermi surface of skyrmions emerges 
and the Landau damping from their particle-hole excita- 
tions gives rise to the dynamical critical exponent z = 3 
for spin fluctuations. As a result, we predict that the 
Gruneisen ratio, given by the ratio between the ther- 
mal expansion and specific heat coefficients, changes from 
r(T) oc T _1 at the antiferromagnetic (Lifshitz) QCP to 
r(T) cx T~ 2 / 3 in the non-Fermi liquid (skyrmion liquid) 
phase. 

II. SKYRMION LIQUID 
A. 0(3) nonlinear a model 

We start from an effective Kondo-Heisenberg lattice 
Hamiltonian, 

H = ^(ffc - ii)c k<J c h v + Jk / ] 8j ■ Si + J Si ■ Sj, 

k i ij 

(1) 



which shows the competition between antiferromagnetic 
ordering (J) from localized spins (Sj) and the formation 
of heavy electrons from the Kondo effect (Jk) of itinerant 
electrons (c^) and local spins. Since we focus on the 
question how the antiferromagnetic order vanishes due to 
the Kondo-coupling effect, the nonlinear a model is our 
good starting point for antiferromagnetism from localized 
spins fill . [l3 | . A careful point in this mapping is that 
itinerant electrons couple to ferromagnetic fluctuations 
of local spins dominantly because the nesting property 
of the Fermi surface does not exist due to their small 
Fermi surface. Their particle-hole excitations give rise 
to the Landau damping dynamics for such ferromagnetic 
fluctuations via the Kondo-coupling effect, given by 



S eff = S NLsM +S Ko ndo, SnlsM = t^- J dr J d d r j i [d T n(r , r)] 2 + [ Vn(r , r)] 2 j + S B , 
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7i£ / 7^( n ( r ' T ) x i 9 Mr,r)}) J^- (n(r',r') x [0 T ,n(r',r')] 



(2) 



SnlsM is the typical nonlinear a model for the 
Heisenberg-type model [l3[ , where the coupling constant 
g r between antiferromagnetic spin waves and the spin- 
wave velocity c r are renormalized by the Kondo-coupling 
term. Sb denotes a single-spin Berry phase term, which 
originates from the path-integral quantization of a spin 
field in the spin-coherent state representation. In the 
present study we will not take into account the role of this 
Berry phase term. SKondo describes an effective interac- 
tion between ferromagnetic fluctuations with A = Jk/ J, 
where the polarization kernel H(q, ifl) = J^\ q \ shows the 
Landau damping dynamics. 

An important remark in this setup is that the Landau 
damping form seems to be robust beyond the present one- 
loop approximation [T^4l6| , where the existence itself of 
the Fermi surface seems to protect the Landau damping 
form. However, we would like to point out that the ex- 



act expression for the polarization kernel still remains as 
an open question. Within this uncertainty, our problem 
is clarified as follows. What is the role of the nonlocal 
effective interaction for ferromagnetic spin fluctuations 
(SKondo) in the antiferromagnetic quantum phase tran- 
sition described by the nonlinear a model (Snlsm)? 
B. Duality transformation 

Performing the duality transformation for this effec- 
tive nonlinear a model, we derive an effective field the- 
ory for skyrmion excitations [17| . First, we rewrite Eq. 
(2) in terms of bosonic spinons given by n(r,r) = 
hzl(r, T)tT aa iz a i (r, t). Second, we consider an easy- 
plane approximation, given by z a — ^e*^", where the 

0(3) symmetry is reduced to 0(2). Resorting to this CP 1 
representation with the easy plane anisotropy, we obtain 
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~ lo ^ I la' ^ I d2r ' CrNFX2 &M r > T ) ~ 9 rM r > t)]U(t ~r',T- r')[d T ^ T (r', r') - d T ^(r', r')],(3) 

I 

where e a is an internal electric charge for bosonic spinons, term gives rise to complications. The resulting skyrmion 
which results from the U(l) gauge redundancy in the CP 1 field theory is given by 
representation, and r is scaled into c r r with f3 r = c r f3. 

It is straightforward to perform the duality transforma- 
tion for Eq. (3) although the Kondo-fiuctuation induced 



S sk = J T dr J dV{ / i sfc $t(9 T -*c T )$ s + |(^ I - lC ^$ s | 2 + m2|$ s |2 + ^|c i , s |4 + _L (axc) 2| 



dr / d 2 r I dr' / d 2 r '— [d x c(r, T)] T U rr ^ TT ,[d x c(r',r')] r , 



(4) 



r 



where $ s (r,T) represents a skyrmion field and Cu(r, r) 



expresses 

2c T N F g 2 r \ 2 



spin-wave excitation. 



and 



J_ 

2 7 <?2 

„ have been done to translate two cou- 
pling constants, J and Jk of the original model into 
two kinds of internal charges of skyrmions, q r and p r . 
Tl rr , TT , represents the effect from Fermi surface fluctu- 
ations, given by H(q,ifl) = |0|/|q| in the momentum- 
frequency space. We would like to emphasize that our 
formulation can be constructed in any general expression 
for the polarization kernel. In this respect we can give 
an answer beyond some limited approximations for the 
polarization kernel. See appendix A for the derivation 
from Eq. (3) to Eq. (4). 

Several remarks should be given for our effective field 
theory [Eq. (4)] of skyrmions. Since two species of 
bosonic spinons exist in the CP 1 representation, there 
must be two kinds of vortices, and identified with 
meron excitations [171 ]. However, we take the limit of 
e a — > oo, resulting in the fact that all gauge non-singlet 
excitations are confined to disappear from the effective 
field theory [ij,[l9j]. As a result, $ s - $ T $| arises nat- 
urally from the confinement ansatz of e a — > oo, identi- 
fied with the skyrmion field instead of meron excitations. 
The effective field theory, Eq. (4) consists of "normal" 
skyrmions and spin wave excitations in the easy plane 
approximation, where the U(l) gauge symmetry is as- 
sociated with the conservation of the z— component of 
localized spins. Another important aspect is the pres- 
ence of the particle-hole symmetry-breaking term, given 
by the linear time-derivative with /z s fe in the skyrmion 
dynamics. As discussed in the introduction, our physical 
picture is based on the analogy with the vortex phase 



diagram of type II superconductors [8(. We speculate 
that there exists a skyrmion "lattice" —like region at least 
at finite temperatures before the antiferromagnetic order 
disappears, which breaks the particle-hole symmetry for 
skyrmion dynamics as the vortex lattice phase. Appear- 
ance of the particle-holc symmetry breaking is attributed 
to the low Neel temperature and the Kondo-coupling ef- 
fect with itinerant electrons 20] . Our physical picture is 
that skyrmions serve effective magnetic fields to itinerant 
electrons, which can play a role in reducing the kinetic 
energy of itinerant electrons. As a result, the total en- 
ergy can be more lowered, where skyrmion fluctuations 
become more softened by the interplay between itinerant 
electrons and skyrmions. Indeed, a similar phenomenon 
has been observed theoretically in the Kondo system on 
geometrically frustrated lattices, where coplanar order- 
ing of localized spins in the absence of itinerant electrons 
becomes unstable to turn into spin chiral ordering struc- 
tures at special fillings of itinerant electrons, which gen- 
erate effective magnetic fields to itinerant electrons and 
quench the kinetic energy to lower the total ground-state 
energy [2lj . The antiferromagnetic QCP is achieved by 
proliferation of skyrmion excitations, but we claim that 
there can exist an intermediate skyrmion liquid state be- 
fore the skyrmion condensed phase. 



An interesting feature of Eq. (4) is the nonlocal term 
with the polarization kernel in the gauge or spin-wave 
propagator. This gives rise to the fact that the "charge" 
p r does not renormalize at the QCP when the linear time- 
derivative term is neglected. It is straightforward to ob- 
tain renormalization group equations for two coupling 
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constants, 



dql 
dhifi 



<?r(e - Vc), 



dpi 
din /i 



P 2 r (e- Vc ). (5) 



? ' c = (TSiLi = ^ c 1r 1S the anomalous dimension of the 
gauge field, where Z c is the wave-function renormaliza- 
tion constant for the gauge field and C c is a positive nu- 
merical constant, e = 3 — d is with dimension d. A key 
point in these renormalization group equations is that the 
anomalous dimension r\ c appears in the same way, which 
originates from the gauge invariance. As a result, we ob- 
tain r\ c = e at the charged fixed point, which is exact 
for all orders [22j . This defines the critical coupling con- 
stant q%. However, there is no equation to fix p r , which 
remains unrenormalizcd. This originates from the nonlo- 
cal interaction generated from the Kondo-coupling effect. 
The non-renormalization for p r motivates us to search 
an appropriate description for the existence of a critical 
line instead of a critical point. One way is to fermion- 
ize skyrmion excitations. The fermionization procedure 
gives rise to non-relativistic dynamics for such fermion- 
ized skyrmion excitations, characterized by z = 2. The 
z = 2 dynamics seems to be inconsistent with the bosonic 
description with 2 = 1. Our requirement of mathemat- 



ical consistency between fermionic and bosonic descrip- 
tions for skyrmions forces us to introduce the linear time- 
derivative term in Eq. (4). The physical picture for this 
particle-hole symmetry breaking term was discussed in 
the previous paragraph. 



C. Fermionized skyrmions and Lifshitz transition 

Possibility of statistical transmutation for vortices has 
been discussed in the vortex liquid phase @ and at quan- 
tum criticality in geometrically frustrated spin systems 
[Io| . Fermionization of skyrmion excitations can be per- 
formed in the same way as that of vortices. The key 
point is that the Chern-Simons term becomes irrelevant 
at quantum criticality 0, [13] , implying that the quantum 
statistics for skyrmions or vortices may not be well de- 
fined at quantum criticality. Actually, such excitations 
are strongly interacting at criticality, where we are not 
allowed to pin down elementary excitations clearly. See 
appendix B for irrelevance of the Chern-Simons term at 
quantum criticality. We reach an effective field theory of 
fermionic skyrmions, 




where ip s represents the fermionic skyrmion field with 
the skyrmion chemical potential ^ and the band mass 
m s k cx fi s k ■ The antiferromagnetic transition is described 
by tuning the chemical potential from p^ < to p^ > 0. 
In this respect the antiferromagnetic QCP is realized at 
(if — 0, which means that the chemical potential touches 
the lowest position of the band, thus identified with the 
Lifshitz transition. The Lifshitz transition is character- 
ized by the dynamical critical exponent z = 2, where 
the renormalization effect from gauge fluctuations is not 
relevant. As a result, both z = 2 critical skyrmion fluc- 
tuations and z = 1 gapless spin fluctuations coexist to 
contribute to thermodynamics. Increasing the skyrmion 
chemical potential from zero to positive, the skyrmion 
Fermi surface will arise. Then, Landau damping for 
gauge fluctuations emerges from particle-hole excitations 
of fermionic skyrmions, resulting in the z = 3 dynamics 
[23l ]. The renormalization group analysis has been per- 
formed [24I [25[ | , where the renormalization group equa- 
tion for q r remains essentially the same as that of Eq. 
(5). An effective field theory has been constructed right 
at the interacting charged fixed point in the absence of p r 



[2a ] , where the coupling constant p r can be easily shown 
to be irrelevant due to the z — 3 gauge dynamics. Irrele- 
vance of p r is consistent with enhancement of the Kondo- 
coupling effect in terms of original variables, implying 
that the skyrmion liquid with a skyrmion Fermi surface 
is a possible solution which originates from the interplay 
between itinerant electrons and skyrmion fluctuations. 

D. Experimental implication: Gruneisen ratio 

Although the skyrmion liquid state can be character- 
ized by the power-law correlation function of skyrmion 
densities, there are more directly measurable thermody- 
namic quantities. Both specific heat and thermal expan- 
sion coefficients can be deduced from two critical expo- 
nents, the dynamical critical exponent z and the correla- 
tion length exponent v. Unfortunately, these two thermo- 
dynamic functions depend on dimensionality. Although 
our skyrmion description was constructed in two dimen- 
sions, we cannot exclude the possible existence of the 
skyrmion liquid state in three dimensions. In this respect 
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it is more reliable to consider a dimension-independent 
quantity, which is the Gruneisen ratio T(T), given by the 
ratio between the thermal expansion and specific heat 
coefficients. It is straightforward to see T(T) oc T~ 1/Vz 
from the scaling argument [26j . 

z v F{T)ocT- 1/uz 

AF-QCP 2 i72 T^ 1 
NFL 3 1/2 r~ 2/3 

TABLE I: Gruneisen ratio at the antiferromagnetic quantum 
critical point (AF-QCP) and in the non-Fermi liquid phase 
(NFL) 

Table I shows how the Gruneisen ratio changes from 
the antiferromagnetic QCP to the skyrmion liquid phase. 
As the dynamical critical exponent changes from z = 2 
at the antiferromagnetic QCP (Lifshitz transition) to 
z = 3 in the non-Fermi liquid phase (skyrmion liquid 
with its Fermi surface), the Gruneisen ratio will follow 
from ~ T _1 to ~ T~ 2 / 3 , respectively. We recall that 
the T -1 behavior was seen from CeNi2Ge2, regarded as 
a conventional antiferromagnet, while the T~ 2 / 3 behav- 
ior was observed from YbRh 2 (Sio.95Geo.o5)2, expected to 
belong to a different class of heavy fermion systems [13] ■ 
In particular, the T -2 / 3 behavior has been attributed to 
the Kondo breakdown mechanism, which differs from the 
present scenario [H, . 



III. DISCUSSION 

A. Role of Berry phase 

The role of Berry phase has been discussed intensively 
in insulating antiferromagnets [l7| . Although it does not 
play any important role in an antiferromagnetic phase, it 
assigns a nontrivial quantum number to skyrmion excita- 
tions in a paramagnetic phase, the condensation of which 
results in translational symmetry breaking (valence bond 
solids). Furthermore, interplay between Berry phase and 
interaction has been proposed to allow exotic excitations 
at a quantum critical point between an antiferromagnetic 
phase and a translational symmetry broken paramagnetic 
state. One dimensional spin-liquid features have been 
claimed to appear in two dimensions, where spin-1 ex- 
citations in both antiferro- and para- magnetic phases 
become fractionalized into spin-1/2 spinon excitations. 
However, the role of Berry phase has not been clarified 
in metallic antiferromagnets. We believe that this prob- 
lem deserves to investigate intensively near future. 



B. On the easy-plane approximation 

To take the easy-plane limit reduces the 0(3) spin sym- 
metry to 0(2) x Z2. Although it is certainly meaningful 



to study this type of antiferromagnets coupled to itiner- 
ant electrons, an actual question is how much different 
physics will arise. 

First of all, we would like to point out that qualita- 
tively distinguished physics between 0(3) and 0(2) sym- 
metries are not kept in the present scheme of approxi- 
mation, where our calculational tool is basically based 
on the scheme of the 1/N expansion. Since critical ex- 
ponents depend on N, one may claim that the quantum 
critical point of N — 2 differs from that of N = 3. How- 
ever, this difference is just quantitative in the 1/N ap- 
proximation. Both cases of N = 2 and N = 3 show 
continuous phase transitions. On the other hand, if one 
claims the emergence of fractionalized excitations in the 
0(N) vector model, possible effective field theories will 
be given in terms of fractionalized excitations (spinons) 
and gauge fluctuations [17] . In this case the symmetry 
reduction can give rise to unexpected errors sometimes, 
for example, resulting in weak first order transitions in- 
stead of expected second order ones [3(J. However, such 
a reduced symmetric model can be manipulated to allow 
qualitatively similar features, compared with the origi- 
nal symmetric model. In particular, deconfined quan- 
tum criticality has been predicted within the symmetry- 
reduced model, which also appears in the full symmetric 
model [l7]. We would like to emphasize that these ex- 
otic excitations are not taken into account in the present 
study. Instead, we consider "conventional" spin 1 ex- 
citations only. In this respect we expect no qualitative 
different physics between fully symmetric and symmetry 
reduced spin models. 



C. More on particle- hole symmetry breaking in 
skyrmion dynamics 

It is our essential proposal to replace the second-time 
derivative term with the linear-time derivative. Unfortu- 
nately, we cannot prove the emergence of the Galileian 
invariance (the linear-time derivative) from the rclativis- 
tic invariance (the second-time derivative) at present. In 
this respect one can call the emergence of the Galileian 
invariance in the presence of itinerant electrons as our 
speculation. 

There must be an underlying physical mechanism 
for this linear-time derivative term, which breaks the 
particle-hole symmetry for "bosonic" skyrmion excita- 
tions. It has been demonstrated that the density of 
skyrmion excitations is finite at an antiferromagnetic 
quantum critical point without itinerant electrons, i.e. 
in an insulating antiferromagnet [311 ] . Of course, only 
the second-time derivative term is allowed in this case, 
i.e., particle- hole symmetric, which means that an equal 
number of skyrmion and anti-skyrmion excitations exists 
at this quantum critical point. Our problem is what hap- 
pens on the particle-hole symmetry in skyrmion excita- 
tions at the quantum critical point where itinerant charge 
carriers are introduced. This is a long-standing problem, 
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where non-perturbative effects from interactions between 
itinerant electrons and "many" topologically nontrivial 
excitations should be taken into account on equal foot- 
ing. Frankly speaking, we do not have any reliable math- 
ematical tools for the description of such interactions. 

Our speculation is that the presence of itinerant elec- 
trons will induce the particle-hole symmetry breaking 
in the skyrmion sector because itinerant electrons favor 
skyrmion excitations (or anti-skyrmions, i.e., one of the 
two). The physical mechanism is as follows. When itin- 
erant electrons move in the background of skyrmions, 
they feel an effective magnetic flux, which quenches the 
kinetic energy of electrons. Our expectation is that the 
gain in the electron kinetic energy contribution is larger 
than the cost in the particle-hole symmetry breaking of 
the skyrmion sector. As a result, the particle-hole sym- 
metry breaking is favorable in the respect of the total en- 
ergy. Actually, this mechanism has been realized in the 
system of frustrated magnets, where the presence of itin- 
erant electrons leads the co-planar ordering in the trian- 
gular antiferromagnet to be ordered into an out-of-plane 
way, which corresponds to a spin chiral order [21| . In 
the present study we do not consider this time reversal 
symmetry breaking because our skyrmion excitations are 
strongly fluctuating at our quantum critical point, which 
will not form such an exotic static order. 



cause we need to calculate spin-spin correlation functions 
of ferromagnetic fluctuations and this is not easy to per- 
form, in particular, for transverse spin fluctuations. The 
spin-correlation function for the z— component is basi- 
cally given by the flux correlator in terms of (V x c) T . 
However, the spin-correlation function for transverse 
components is related with "magnetic monopole" corre- 
lations [10| which changes 27r-flux in (V x c) T because 
transverse components flip the z— component spin and 
the z— component spin is identified with the gauge flux 
in the dual representation. 



F. Comparison with a non- Fermi liquid metallic 
phase in the slave-fermion theory 

We would like to point out that a possible non-Fermi 
liquid metallic state has been prop osed in the context 
of the slave-fermion theory [29|, l32j |. However, the non- 
Fermi liquid state from the slave-fermion theory is char- 
acterized by the existence of spin gap while the skyrmion 
liquid state does not exhibit such a spin gap, where spin 
fluctuations remain gapless. The skyrmion liquid phase 
should be distinguished from the spin-gapped non-Fermi 
liquid state. 



D. Summary of approximations 

It is necessary to summarize our approximations for 
the procedure from Eq. (1) to Eq. (6). It is basically 
an exact procedure to derive an effective 0(3) nonlin- 
ear (7 model [Eq. (2)] from the Kondo-Heisenberg lattice 
Hamiltonian [Eq. (1)] after integrating over itinerant- 
electron degrees of freedom if we replace the Landau 
damping form for the polarization kernel with its gen- 
eral expression. Although we keep the Berry phase term 
in Eq. (2), we do not take into account its role in spin dy- 
namics, particulary, skyrmion dynamics. The easy-plane 
approximation has been resorted for the derivation from 
Eq. (2) to Eq. (3). The effective dual Lagrangian [Eq. 
(4)] in terms of bosonic skyrmions has been found from 
the duality transformation, as shown in appendix A. An 
essential aspect is to replace the second-time derivative 
term with the linear-time derivative. This is our key 
proposal. The resulting field theory [Eq. (6)] in terms 
of fermionic skyrmion excitations come from the bosonic 
skyrmion theory [Eq. (4)], performing the fermionization 
procedure. In appendix B, we prove irrelevance of the 
Chern-Simons term at quantum criticality, which sup- 
ports the emergence of fermionic skyrmions. 



E. On electron dynamics 

It is not straightforward to understand the feedback ef- 
fect from skyrmion excitations to itinerant electrons be- 



IV. CONCLUSION 

In conclusion, we proposed a quantum critical metal- 
lic state for a non-Fermi liquid phase of Ge-substituted 
YbRli2Si2, identified with a skyrmion liquid state 
with a skyrmion Fermi surface. Nonlocal interactions, 
originated from Kondo fluctuations, and fermionized 
skyrmions, argued to be allowed from both nonlocal in- 
teractions and coexistence with itinerant electrons [9l[l0|. 
give rise to a stable interacting fixed line. The antifcrro- 
magnetic transition has been interpreted by the Lifshitz 
transition of fermionized skyrmions, characterized by the 
dynamical exponent z = 2. The skyrmion liquid state is 
described by its Fermi surface, resulting in z — 3 for col- 
lective spin fluctuations. We suggested the fingerprint 
of this scenario from the change of the Gruneisen ratio, 
where it behaves from ~ T -1 at the antiferromagnetic 
quantum critical point (z = 2) to ~ T~ 2 / 3 in the non- 
Fermi liquid phase (z = 3). 

An important aspect of the present study is on the 
introduction of Fermi-surface fluctuations. Generally 
speaking, it is not justified to integrate over Fermi-surfacc 
fluctuations. An interesting point is that although we 
cannot estimate the precise form of the polarization ker- 
nel, which results from Fermi-surface fluctuations to give 
nonlocal interactions between collective spin fluctuations, 
such nonlocal interactions can be irrelevant if we assume 
that skyrmion excitations become fermions instead of 
bosons with particle-hole symmetry breaking. 

The next question will be, "Can we derive this physics 
based on a purely diagrammatic way?" The present re- 
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suit can be translated into the fact that nonlocal inter- 
actions from Fermi-surface fluctuations may cause inter- 
esting novel physics to dynamics of spin fluctuations. Al- 
though this issue has been addressed in the duality pic- 
ture, it is an important direction to prove the emergence 
of z = 3 antifcrromagnetic quantum criticality based on 
the original description for spins. 
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Appendix A: Duality transformation 



Starting from Eq. (3) and performing the Hubbard-Stratonovich transformation, we obtain 



7 = j DMr, 
-ij T (8 T <p t - r 04.)} 



T)Da^{r, r)DJ°(r, r)Dj T (r, r) cxp 



dr 



J? -iJ£(0 M 6x-a„) + ^2(0xa) 



(It j (f-r \ &T dr 1 I d 2 r 



o 



1 



4c r N F X 2 



j T {r,T)[Tl{r - r' ,t - T')]-^j T {r' ,t') 



(Al) 



where J° and j T are auxiliary fields which correspond to currents. 

Separating the phase field into spin-wave (smooth) and vortex (singular) parts and integrating over spin-wave 
fluctuations, we obtain equations of constraints, 



dp, 31 + vd^juS^ = 0, 
where a = ± represent spin f and I, respectively. These equations are solved to give 

where is U(l) gauge field to incorporate spin-wave excitations. 



(A2) 
(A3) 



Integrating over J° with Eq. (A3), we obtain 
Z = J Dcf> v cr (r,T)Da^r,T)Dj r (r,T)Dc°(r,T)exp[- Jj dr j d 2 r[^[(8 x - CT j M V] 2 - W * - <v) 



(A4) 



It is straightforward to perform the integration for j T . Then, we obtain a dual action in the first-quantization 
expression 

Z = j D4>l(r,r)Da^r,T)Dc^(r,T)e^[- j\r j d 2 r{^-(8 x c°f - i{8 x c<%(c^ - a M ) + ^(8 x a) 2 } 
~ So ^ J Jo r ^ ' J d2r '^[^ xct ^ _ ( axci ^]( r ' T ) 

{g r 6(r - r')8{r r 1 ) + [U(r -r',r- r'T^Kd x c^) T - (8 x c+) r ](r', /)] . (A5) 

We note that dynamics of spin-wave excitations is strongly modified by nonlocal interactions from Fermi-surface 
fluctuations of itinerant electrons through the Kondo-coupling effect. 

Considering the minimal coupling term of c°{J v ]° in Eq. (A5), where {J v ]° = (8 x dcf)^)^ is identified with a vortex 
current, one can construct the second-quantization form of Eq. (A5), introducing vortex field variables. As a result, 
we obtain an effective dual field theory 

2 + ^l<f.| 2 + yl^| 4 + f(3x0 2 



Z = J D<t> a (r,T)D afl {r,T)Dc°(r,T)ex.p[- J " dr J d 2 r(\(d^ - ic°)<£ a 
+i(8 x + ±(d x a) 2 ) - dr J d 2 r dr 1 j d 2 r / c r N F g 2 X 2 [(8 x c^) T - (8 x c+) r ](r, 



II(r -r',T- t')[(8 x c*) T - (8 x c^) T ](r',r') 



(A6) 
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$o- represents a meron field, which acts as a source of (spin- wave field), where the meron current is given by 
[J v ]^ = —i[&[ r (d ll $ (T ) — (9 M $J.)$o-]. m 2 is a mass of a meron, which identifies its quantum critical point with m 2 = 0. 
u v is a coupling constant for its local self-interaction, determined from non-universal short-distance physics. 

As discussed in section II-B, we will not allow meron excitations. Emergence of meron excitations has been discussed 
intensively in Ref. TtJ , where the interplay between Berry phase and local self- interactions of merons would make 
"magnetic" monopole excitations of compact U(l) gauge fields a, become irrelevant in the renormalization group 
sense, giving rise to merons as elementary excitations. Emergence of meron excitations implies that of spinons, which 
carry fractional spin quantum number 1 /2. On the other hand, we take the limit of e a — > oo, where such fractionalized 
excitations cannot exist. Merons should be combined with anti-merons, forming skyrmions, which means that spinons 
should be confined with anti-spinons, resulting in spin 1 excitations, identified with conventional spin fluctuations. It 
is still extremely difficult to clarify the precise condition for the emergence of deconfinement although its existence 
seems to be accepted. 

In order to realize the confinement ansatz, i.e., e a — > oo, we take the amplitude-frozen limit of Q a ~ e l6 " and 
obtain 

(W* ~ <£ f = \ ([dA + d^] [ C t + <#) 2 + \ {[d^ d^] [ C t - eft) 2 . (A7) 
At the same time, we rewrite the kinetic energy of the vortex gauge field (spin-wave) as follows, 

|(S x c°f = |(p x ct) + (d x c^)]) 2 + |(p x c t) - (d x ct)])'. (A8) 
Then, one can rewrite Eq. (A6) in terms of newly defined vortices and gauge fields, 

C v = \(d, zc+)$+| 2 + m 2 + |$ + | 2 + ^t|$ + | 4 + ia^ vX d v c+ + |(d x c+) 2 + i(9x a) 2 



-|(5, - lC -)<£^ + ml\*>_\ 2 + ^|$_| 4 + |(9 x c-f 

fir 



+c r N F g 2 X 2 ' dr' J d 2 r'[dxc-{r,T)] T Ii{r-r',T-T')[dx.c-{r',T')] T , (A9) 

where 9± — — > $± is a "new" vortex field and = - JL ^ iL is also a "new" vortex-gauge field. An essential 
aspect of Eq. (A9) is that only the $+ sector contains interaction with compact U(l) gauge fluctuations a, while the 
dynamics of $_ decouples from the gauge dynamics completely. If we take e a — > oo, $ + cannot appear as elementary 
excitations. This is the confinement ansatz, which do not allow such fields as carry the U(l) internal gauge charge 
of a M . This dual language can be interpreted as the fact that either spinons or anti-spinons cannot appear as the 
physical spectrum. 

Integrating over a, in the limit of e a — > oo, we see that $+ and c+ disappear in the physical spectrum. As a 
result, we reach an effective field theory for skyrmion dynamics in the presence of itinerant electrons 



c sk = \(d, - lc ;)$_| 2 + m 2 ^! 2 + ^rl<f-l 4 + f (3 x c- f 



+c r N F g 2 X 2 J J dr' J d 2 r'[d x c"(r, r)] T U(r - r', r - T ')[d x c"(r', t')] t . 

Appendix B: Irrelevance of the Chern-Simons term at quantum criticality 

Starting from Eq. (4) and replacing $ s with 1 & s , we obtain 



(A10) 



Z sk = [ D9 a Dc^exp\- ^ dr [ d 2 rU\ (d T - ic r )<f> s + — - ic^ s \ 2 + ^\<t> s \ 2 + ^-|1> s | 4 + ^(d x cf 

J L Jo J L Ms/c Msfc ^^sk Z 1r 



J2 / 



dr I d 2 r I dr' I d 2 r'^[d x c(r,r)] T Il rr , iTT ,[fl x c(r',r')] T ■ (B 
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2/4 



In order to introduce the statistical transmutation for skyrmion excitations, we perform the Chern-Simons gauge 
transformation, which attaches a flux to a skyrmion. Then, we obtain 

Z sk = ( D^sDc^Da^cxp - I dr f d 2 ri^l(d T -ic T -ia T )^ s + —\(d^-ic^-ia^ s \ 2 + ^\^ s 

+ 7^(d x cf - i-L a ^ x d„a x } - J\t J d 2 r J\t' J A^[axc(f,r)] f n rr , rr -[axc(r',/)] r ], (B2) 

where ^ s is a fermionized skyrmion field and is the Chern-Simons gauge field to guarantee the statistical trans- 
mutation. 8 is the Chern-Simons angle, where 9 = ir describes the fermion-boson or boson-fermion transmutation. 
Shifting the skyrmion gauge field (spin wave) as c p — > c p — a M , we obtain 



Z sk = I D^sDc^Dav expf- dr [ d 2 r\¥ s (d T - ic T )V s + — - ic^ s \ 2 + 

J 1 Jo J 1 Msfc ^sk 

I i 

+ 2q 2 ^ X ° ~ 9 X ~ i ^ a ^ e ^ xd ^ ax \ 



(B3) 



- J dr J d 2 r dr' J d 2 r' ^ [9 x c(r, r)-3xa(r, T)] T IW, rr , [9 x c(r', r') - d x a(r', t% 
Integrating over the Chern-Simons gauge field, we reach the following expression 

f" dr f d 2 r{¥ s (d T *c T )* s + -j-\(d, - ic,^ s \ 2 + ^|* s | 2 + ^|* s | 4 + ^ (d x c) 2 } 

JO J *■ Ms/c Msfe Z / i sfc Z< ?r 



L»* s L>c M exp 



0r 



rfr y rf 2 r y dr'y d V^[9 x c(r,T)} T Il rr ,, TT ,[d x c(r',T%-S irr .[ Cft {r,T) 

where Si rr . [c M (r, r)] describes the dynamics of c M , which results from the integration of the Chern-Simons gauge field. 
It turns out to be 



(B4) 



<5irr.[c M (r,r)] oc c fl (q,ifl)U^(q, ito)c v (-q, -ifl) cx \q\ s c^{q,iQ)P^ v {q,iQ)c v {-q, -ifi), 



(B5) 



where P^ v {q, ifl) is the projection operator to the transverse direction. As shown in this expression, the proportionality 
to |q| 3 implies the irrelevance of the Chern-Simons interaction at quantum criticality. 
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